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Abstract 



We explore the possibility of detecting flavor changing neutral Higgs couplings 
at the Next Linear Collider (NLC) through e^e~ — > V(.i^etc. In the framework 
of a general two-Higgs doublet model, we perform a complete calculation and 
find that a (e^e~ — > v^Vgic, Ve^^eic) could reach ~ 9 fb for ^/s = 2 TeV. This 
amounts to an annual production of 500 tc plus tc pairs at the NLC with 
an integrated luminosity of 50 fb~^. The dependence of tc-production rate 
on the neutral scalar mixing angle is mild except when sin^ a — > or 1. The 
uuW^W~ background should be manageable after 6-tagging, while vuti back- 
ground should not be a problem when the signal event rate is still interesting. 
The process, together with e+e~ — > i>ei'eW~^W~ , V(,V(^ZZ studies, offer the 
chance of measuring the t-c-Higgs coupling. 

PACS numbers: 14.80. Cp, 14.65.Ha, 12.15.Ff, 13.90. +i 
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I. INTRODUCTION 



The mechanism for symmetry breaking and the fermion mass and mixing hierarchy 
pattern are the two remaining mysteries in the electroweak theory. The construction of high 
energy colhders such as the LHC and NLC are in fact aimed at resolving such mysteries. In 
this regard, the physical processes that should be studied thoroughly at such machines are 
those involving the top quark, whose properties have yet to be studied carefully, as well as 
the yet to be discovered Higgs boson(s). 

It was suggested some time ago that [||, in multi-Higgs doublet models, the "natural" 
flavor conservation condition is not mandatory for the suppression of flavor changing 
neutral current (FCNC) processes. Rather, Nature has provided its own cure: the ex- 
isting hierarchical patterns in quark masses and mixing angles may imply a pattern for 
flavor changing neutral Higgs couplings (FCNH) that is consistent with low energy data 
||T|]. An interesting consequence of this framework is the possibility of sizable t-c- neutral 
Higgs couplings which would have notable impact on top quark and Higgs physics 0,0. 
To probe such couplings at colliders, several processes have been proposed which can 
be studied at the NLC or LHC. At the NLC, one may look for tc pair production via 
e+e" — Z* — >• tc, ic [|] (where the Z-t-c coupling is loop- induced), or like-sign top pair pro- 
duction via e+e" h^A^ ttcc, ttcc [§. At the LHC, such flavor non-diagonal couplings 
can be probed through the parton subprocess eg — >• tA'^ — > ttc [0 , which involves the FCNH 
coupling directly in the production process. 

Recently, Bar-Shalom et al. pointed out that FCNH couplings may be probed at 
the NLC via the WW fusion process e^e'^ — > Uei^etc, UgUetc, as shown in Fig. 1. With 
\/s = 2 TeV, and the masses of neutral Higgs bosons being 250 GeV and 1 TeV respectively, 
they found a^utc = (e"'"e~ Uei'etc) + a {e~^e~ VeT^etc) 5 fb. Assuming an integrated 
luminosity of 50 fb~^ at the NLC, this implies an annual production of 125 tc and an equal 
number of ic pairs. The process has a much larger tc production rate than e^e^ Z* ^ tc, 
and does not suffer from s-channel suppression as e^e~ h^A^ ttcc. In view of this, we 
would like to follow up on this work. We shall perform a full calculation and compare with 
the effective W approximation used in Ref . [8] , explore different scenarios for neutral Higgs 
masses, and clarify parameter dependence of the tc production cross section. 

The paper is organized as follows: In Section II, we briefly review the two-Higgs doublet 
model with FCNH couplings and present the result of a full calculation of a^utc using helicity 
methods. We then point out that av^tc is largest when both neutral scalars have mass of 
order the weak scale. This becomes the focus of our discussion throughout the paper. In 
Section HI, we demonstrate the utility of the narrow width approximation. In Section IV we 
show that a^ptc is not sensitive to the mixing angle of neutral scalars, and remains at the fb 
level for y/s > 1 TeV. Some discussion of signal vs. background is given. After concluding 
in Section V, we leave some technical details in Appendices A and B. 

II. FULL CALCULATION 

The calculation of ai,utc is based on the Lagrangian of a general two-Higgs doublet model 
with flavor changing neutral Higgs (FCNH) couplings 
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- (uLM^u^ + dLM'^dn) 

+ {-UL&'>un + dL^^''UR) tV2lm<Pl 

- dLV^^^''^URV2(P^ + ULV^^''URV2(Pt +H.C., (1) 

where ul,{r) and dL,{R) are flavor multiplets of up-type and down-type quarks respectively, 
and M^''^ are their diagonalized mass matrices. Note that we have relegated all the FCNH 
couplings to the second doublet as a result of rotating to the specific basis (^g) = and 
(05) = "w/v^ [0. This is because there is no discrete symmetry [|] as in usual two Higgs 
doublet models ||I0[ to distinguish between $1 and $2, so the familiar tan/5 = Vi/v2 pa- 
rameter is not physical. Assuming CP invariance in the Higgs sector, the scalar fields 
v^Im02 and (pf are identified as the physical pseudoscalar boson and charged scalar 
if^. The CP even neutral scalars \/2'Re(f)'( and a/2 Re 02 can still mix through the Higgs 
potential V^($i, $2) into the physical states and h^. In the limit that the mixing angle 
sin a — i> 0, ^ -\/2Re05 becomes the "standard" Higgs boson with diagonal couplings, 
while ^ \p2 Re 0° has FCNH couplings characterized by the non-diagonal matrix . 
In our calculation as well as in Ref. [8], the simple Cheng-Sher ansatz is adopted: 



= h (2) 

where fi/s are constants of order unity. The coupling is expected to be the largest and has 
the most prominent signature to be searched for in collider experiments. From Eq. (1), we 
can now single out the relevant couplings for computing the process e"'"e~ i^e^^etc, Ue^etc 
given in Fig. 1. Since we wish to compare with Ref. [8], we take ftc — v^. The resulting 
t-c-Higgs and Higgs-PF-iy couplings are 

Cint. = ^ (sina -|- cosa /i) tc -|- h.c. + gmw (cosa H — sin ah) W^Wfj^. (3) 

V2mw 

The Higgs-Z-Z couplings can be easily incorporated, and the cross sections for e'^e~ —>■ 
e~^e~tc, e^e~ic via ZZ fusion are simply related to that of e~^e~ — >■ UeUetc, Pei^eic [El. 



A. Helicity Amplitude Calculation 



A full calculation of a^ytc is rather involved as the process considered is a 2 — > 4 scattering. 
An efficient way of doing it is by employing the helicity method [|l^, which facilitates the 
numerical manipulations of Feynman amplitudes. 

The amplitude for e+(pi)e~(p2) ^e{pz)^e{pA)i{pt)c{pc) reads 



iM = F Ai)7^P_i;(p3, A3)] [u{p4, X4)-f^'P_u{p2, X2)] [u{pt, Xt)v{pc, X^)] 



g2 - m'jj + imn'^H q'^ - ml + irnhVy, 



(4) 
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where q is the momentum of the intermediate Higgs boson, P± = (1 ± 75)/2, and 



F = cos a sin a ^-^^ {igmw) ( ^^^^^^" ^ ' (^) 

is a collection of coupling coefficients. Note that, except for the relative sign and differences 
in mass and width, the h and H contributions are basically the same. All fermion masses 
are set to zero except for the top quark, and the rric dependence is kept only in the coupling 
of Eq. (2). The helicities of leptons are therefore completely fixed by their left-handed 
vector couplings to W bosons, i.e. Ai = A3 = + and A2 = A4 = — . However, there are four 
combinations involving the helicities of top and charm quarks. 

Let v4(Ai,A3) = Ai)7^P_w(p3, A3) and ^(Aa, A4) = u{p^, \^)'j''P_u{p2, \2)- One 
finds (see Appendix A for details) 



A{++) = J2E,J2E3 {P3 + |7m-Ipi+) 



B{—) = pE,pE, (p4 - |7^|P2-) , (6) 

where uj±t = \J Et ± 7± are 2x2 matrices defined by 7^^ = (1,=fc?), and |p±) denote 
the two-component eigenvectors of the helicity operator p ■ that is 

, / cosf \ , /-e-^"^sin|\ 

l^+)=(e^^sin|)' l^-^=( cosf 0' 

where 9 and are angles for p. For the top-charm scalar density, define C(Ai, Ac) = 
At)f (pc, Ac), one gets four combinations (see Appendix A for details) 



C{TT) = -pE, uj+t {Pt T \Pc±) 



C{t±) = -pE, u^t {pt T \PcT) . (8) 

Since A{++) and B{ ) are already fixed, there are four helicity amplitudes iM.{Xt, Ac) oc 

C(Af,Ac). With all four helicity amplitudes constructed, the subsequent numerical calcula- 



tions can be done in a straightforward manner by utilizing the program ONETOP |T2 . 



B. Comparison with Bar-Shalom et al. 

To compare with Ref. [8], we compute CFyytc for uih = 1 TeV and sin^a = 1/2. The 
cross section a^i^tc as a function of rrih for ^/s = 0.5, 1, 1.5 and 2 TeV are shown in Fig. 2. 
It peaks notably at rrih ~ 250 GeV and decreases monotonously as rrih increases from 250 
GeV. In accordance with the difference in the propagators given in Eq. (4), it vanishes in 
the degenerate limit nih = mu = 1 TeV. This is a special case for the choice of sin^ a = 1/2 
(i.e. a = vr/4), because the Higgs properties are identical in the degeneracy limit, so the 
amplitudes arising from each Higgs would then cancel completely. For y/s = 2 TeV, the 
maximal value of cri,utc is around 4.5 fb, which is smaller than 5.2 fb obtained in Ref. [8] 
which uses the effective W approximation. Such an overestimation by the effective W 
approximation is a typical phenomenon in collider physics [0. 
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The prominent peaks in Fig. 2 suggest that the cross section arising from h alone 
would be the largest at rrih ~ 250 GeV. Similar behavior should then be expected for the 
contribution from H. We therefore expect the total cross section resulting from H and h 
to be the largest if both niu and rrih are of order the weak scale. This precisely fits the 
arguments given in Refs. [6,7] which emphasizes the mass range 

200 GeV < nih, mu < 2mt ^ 350 GeV. (9) 

The lower bound is to allow the tc threshold to turn on. The upper bound of 2mt ^ 350 GeV 
was imposed originally for the pseudoscalar A^. For h and H, as can be seen from Fig. 2, 
the cross section is still sizable up to ruh, h — 500 GeV for y/s > 1 TeV. This is because 
T{h, H — >• VV) ^ r(/i, H — >• /'/) for ruh, h ~ 350 GeV, and the opening of ti mode does 
not increase substantially the total width of h or H. However, for the range of Eq. (9), the 
ti background to the tc or ic modes would be suppressed. 

To show that a,yi,tc is indeed more significant in the the range of Eq. (9), we show in Fig. 
3 the cross section a^^tc as a function of nih for mn = 300 GeV and sin^ a = 1/2. The cross 
section drops to zero at the degenerate limit rrih = rrtH = 300 GeV in a much more dramatic 
way. However, such a severe cancellation does not generally occur since there is no reason 
for nih and mn to be degenerate, and the cancellation is anyway incomplete for other values 
of sin a. The cancellation effect is negligible if the mass difference AM = \mH — fahl is a few 
times the widths of both Higgs bosons (see Appendix B). Slightly away from the degeneracy 
limit, the cross section rises to its peak value = 8.0 fb at rUh ~ 250 GeV for a/s = 2 TeV, 
which is almost twice as large as the case with itih = 1 TeV. As m/j increases to 1 TeV, 
avvtc drops to about 3.6 fb, which is mostly from H . For a lighter /i, i.e. nih < 250 GeV, the 
cross section also drops. This once again illustrates the fact that cr^jytc receives the largest 
individual contributions from h and H respectively at rrih, mu ~ 250 GeV. 



III. THE NARROW WIDTH APPROXIMATION 

It is important to note that the widths of neutral Higgs in the mass range of Eq. (9), even 
up to ~ 500 GeV, are still quite small compared to their masses. The Standard Model (SM) 
Higgs boson ifsM provides an upper bound to H and h widths, for example Tugyi ~ 15 GeV 



for mHsM ~ 350 GeV [|IOl- Since the widths of both Higgs are narrow in the mass range 
of interest, it is convenient to compute cr^i^tc in the narrow width approximation. We may 
approximate cr,^^tc by the cross section of Higgs production cr(e+e~ — > u^u^h^H)) multiplied 
by the branching ratio of the flavor changing decay h{H) tc, ic. This approach is much 
simpler than the previous full calculation or even the effective W approximation. One can 
then determine the Higgs mass and sin^ a dependences of a^utc with ease. 



A. WW -^h, H Production 

Compared to the previous calculation of a{e^e^ uutc), it is considerably simpler 
to compute the cross section a{e~^e~ — > UeJ^eh{H)). It is identical to that of SM Higgs 
production a{e~^e~ — >• Uei^eHsu) = o-uuHsm [II3l' except for the additional factors of cos^ a or 



sin^ a. The amplitude for e'^e z/gZ/e-f^sM is 
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8 \2pi-qi + J \2p2 ■ q2 + ) 
X [^^(Pi,-5i)7''(l -75)t;(gi,S2)] [M(g2, 54)7/^(1 -75)m(P2, S3)]. (10) 

Averaging over the initial and summing over the final state spins give 

-j:\M? = g'ml. [ L ^2 )^(^2 ■ gi)(Pi ■ g2), (n) 

where we have neglected fermion masses. The final state phase space integration is done by 
VEGAS For = 2 TeV and m^sM = 250 GeV, we find a(e+e- ^ Uet^eHsu) ~ 264 fb. 
The cross section for other values of rriH and ^/s can be read off from Fig. 4. 

'R. h, H ^ tc Decay 

To compute the branching ratio BR(/i, H ^ tc), we note that the dominant decay modes 
for nih^ H < 2mj are h, H ^ W~^W~ , ZZ, bb []TD[ and tc, tc where the latter are specific 



to the current model. The width of each decay mode is well known: 

2 3 

T{H ^ W^W-) = cos^ a - Ax^, (1 - Ax^, + 12x^), 

2 3 

9 m'lj 



T{H ^ ZZ) = ^ \ cos' a^l- Axl (1 - 4x| + 124), 



r(H ^tc + tc) = X m^m, sin^ a (1 - x^)^/^ (1 - x'_f'\ (12) 

with Xi = rrii/mH and x± = {mt±mc)/mH- For T{h W~^W~ , ZZ) , etc. one simply 
changes sin^ a — > cos^ a. Note that we have assumed SM couplings for bb, although it should 
depend on more parameters (this is another reason for the mass range of Eq. (9) so we avoid 
uncertainties in H{h)-t-t coupling). However, the bb mode is unimportant for our purpose. 

For a generic mixing angle a, vector boson decay modes dominate over the fermionic ones 
since the former is proportional to while the latter only depends on rrih^ h linearly. 

One can clearly see in Fig. 1 of Ref. P] this severe suppression of BR(/i, H tc) for generic 
a values [jl4[. However, for extreme values of a ^ or 1, the WW, ZZ modes could be 
very suppressed, and either BR(/i — tc) or BR(if tc) become significant p. 

The threshold behavior of the tc mode and the dominance of h, H WW, ZZ modes 
in general help us understand the peak in a^,ytc at rrih^ h ~ 250 GeV. We show in Fig. 5 the 
mass dependence of BIi{h tc) for a few values of sin^ a in the range 

0.1 < sin^a < 0.9. (13) 

BR(iJ tc) can be simply obtained by making the change sin^ a cos^ a. We do not 
include extreme cases of sin^ a ^ or 1 since cr^j^tc in these limits. The shape of 
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Fig. 5 can be understood as follows. BR(/i, H tc) rises sharply right after the opening 
of the tc production threshold. The growth of BR(/i, H tc) should however stop at 
certain Higgs mass, since T{h, H ^ VV) is in general dominant and grows more rapidly 
as mh,H increases. The peak position = 260 GeV for BR(/i, H ^ tc), which is the 

main reason behind the peaks seen in Figs. 2 and 3, marks the point where the growth in 
T{h, H tc) is overtaken by T{h, H VV). It is interesting to note that BIi{h, H — > tc) 
always peaks at nih, h — 260 GeV independent of the sin^ a we choose. This is easily 
understood since, for generic a, BR{H tc) = T{H tc)/TH ^ T{H tc)/T{H VV), 
i.e. BR(if — > tc) ~ tan^ a fimn) where firnu) is largely a-independent and peaks at 
mu = 260 GeV. Similarly, we have BR(/;, — > tc) ~ cot^ a f{mh). Such a simple dependence 
on the mixing angle a makes Fig. 5 very useful. For any sin^ a in the range Eq. (13), one 
can obtain the branching ratio BR(/i, H tc) for any Higgs mass by simply scaling via the 
relation BR{h, H ^ tc) ^ cot^ a f{m\), tan^ a f{mjj). 

We note that the kink due to ti threshold becomes more visible for small sin^ a values. 
This is because the VV contribution to the Higgs width becomes suppressed and the relative 
weight of the ti contribution becomes more significant [|T^ . Such a kink is not apparent in 
Figs. 2 and 3 because the sin^ a = 1/2 case was used. 

C. Cross Section 

The SM Higgs width provides an upper bound to Th and Th- We can therefore use the 
narrow width approximation for rnn, h < 500 GeV. The cross section auutc can be written as 

(Tuutc = (TuuHsmi'^HsM = ^h) X COS^ « X BR(if tc) 

+ (^vvHsmi'^HsM — ^h) X siu^ a X BR(/i tc) + interference terms, (14) 

where i^sM denotes SM Higgs. Note that, with |m|^ — f^W > (3 — 4) x rrihTh, the interfernce 
term can be safely neglected (see Appendix B for details). 

To locate the peak of a^utc for generic sin^ a, Eq. (14) can be rewritten as 

CTuutc = CTuuHsmi^H) X siu^ a X finin) + CFyyHs^{mh) x cos^ a x /(m/^), (15) 

where we have neglected the interference term by assuming a large enough splitting between 
rriH and rrih. With niH and sin^ a fixed as in the case of Figs. 2 and 3, ay^tc only depends on 
o'uuHsMi''^h) X firrih). Since fijrih) peaks at rrih = 260 GeV and cTuuHsm is a monotonously 
decreasing function of itihsm^ position of nih giving maximal auutc should be shifted 
downward from 260 GeV. This is exactly the case as seen in Figs. 2 and 3 where such effect 
are most significant for y/s = 0.5 TeV since ay^Hsm drops most steeply for increasing m^^g^ 
for this case. For ^/s = 2 TeV, this shift becomes much smaller as (TwHsm is relatively flat. 

IV. DISCUSSION 

To illustrate our arguments so far, let us explore the "maximal" and "minimal" (Jpjjtc cross 
sections in the mass range of Eq. (9), and sin a and -/s dependences. We shall also make 
some general discussions about signal vs. background and compare with other processes. 
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A. Range of Cross Sections 



For "maximal" a^^tc, take , for example, uih = 250 GeV and rrih = 240 GeV so — 
mil > 4 X rrihTh, and the interference term in a^^c can be safely neglected. Since the masses 
are approximately equal, one can rewrite Eq. (14) as 

Ci^i/tc ~ (^uvHsM X (cos^ a BR(i7 tc) + sin^ a BR(/i tc)^ , (16) 

where the mass of Hsu can be taken as either that of H or h. Note that the combination 
cos^ a BK{H tc) + sin^ a BIi{h tc) determines the sin^ a dependence of (Tuutc, which is 
plotted in Fig. 6. It is interesting to see that both cos^ a BK{H tc) and sin^ a BR(/i tc) 
are sensitive to sin^ a but their sum is not. This is in large part because we chose almost equal 
rriH and m^, and reflects the mutually compensating nature between the two contributions. 
The effective fraction cos^ a BR{H tc)+sm'^ a BR{h tc) of the ("SM") Higgs production 
cross section stays between 2 — 3% for almost the entire range of sin^a of Eq. (13), but 
becomes extremely suppressed for sin^ a outside this range.. 

For = 2 TeV, sin^ a = 1/2, and tuh, ruh = 250, 240 GeV, from Figs. 4 and 6 we find 

(T^^tc ~ 270 fb X 3.2% = 8.6 fb. (17) 

This is in good agreement with the maximal cross section obtained earlier from the full 
calculation, and illustrates the effectiveness of the narrow width approximation. The sin^ a 
dependence is very mild. For example, at sin^ a = 0.1 or 0.9, a,^i,tc — 6.8 fb for ^/s — 2 TeV, 
which is still comparable to the maximal cross section. The sin^ a dependence for individual 
h or H contributions is much more significant. 

To explore the "minimal" cross sections within the range of Eq. (9), we note from Figs. 
2 and 3 that the contribution of h, H to Cuytc is roughly equal for nih, h = 200 GeV and 
nT'h, H — 350 GeV. We therefore present the results for uih — 350 GeV and nih — 200 GeV, 
which gives roughly the smallest (Ti,i,tc ior the mass range of interest. Wc plot in Fig. 7 
Cuvtc for this set of Higgs masses as a fuction of sin^a for ^/s = 0.5, 1, 1.5, 2.0 TeV. It 
is seen that sin^a dependence remains mild. What is remarkable is that, for almost all 
values o/sin^ a, a^^tc is at fb level or higher for ^/s > 1 TeV. This promising result for cr^j/tc 
holds only in the mass range given by Eq. (9), although the range can be extended up to 
mh, H ~ 400-450 GeV or so. 

In both Figs. 6 and 7 we have illustrated with cases where the h and H peak (in sin^ a) 
contributions are comparable, hence their sin^ a dependences are mutually compensating. 
For more general choices of riih and trh values, some sin^ a dependence would remain for 
(^vvtc-i which is reflected in and easily scaled from the individual h or H contributions. 

B. Signal vs. Background 

Turning to the experimental signal at the NLC, one needs to consider the final states 
from top decay, t — > 6jij2, hence the signal modes are 

e^e~ > Ve^etC > Ve^e^i i'^bc, V^Ve bcjij^, (18) 
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and similarly for v^Vi^c. Since typical cross sections are a few fb in the mass range of Eq. 
(9), with an integrated luminosity of 50 fb~^, we expect of order 100 or so (no more than 
300) vvhcj\22 events, and 1/6 of this in each + Ihc channels, where £ = e^, /z^, r^. 
Although the event rates are significant, we find that the latter is not very promising once 
backgrounds are taken into account. 

What are the potential backgrounds? Since the z/gZ/g pair should carry away missing 
transverse energy ~ mw-, WW fusion events should be relatively distinct at the NLC 
For the mass range m^^u < 2mt of Eq. (9), background from e~^e~ — > i^e^eh, i^e^^eH —>■ 
I'ei^eti is absent. The major background to be considered is therefore e"^e~ — > Uei^eW^W' 
since it is more abundantly produced via h, H W^W" . From Figs. 4-7 one sees that 
(j(e^e~ Vf-VeW^W') is typically 20 to 30 times larger than the signal. But the important 
point to notice is that W decays do not contain h quarks (< 10~^ in BR). The chief tool to 
suppress the WW background is therefore 6-tagging, expected to be very efficient at Linear 
Colliders \T^\- However, since 1/3 of ly decays contain charm quarks, fake rate of 6-tagging 



might be an issue. In particular, the 'iv + 1 + hc mode would not be easy to distinguish from 
3z/ + £ + cs fakes when the signal event rate is so low. In contrast, the vv bcjij2 mode has a 
second handle: kinematics and full reconstruction. With one 6-tagged jet, two of the three 
remaining jets should reconstruct to mw and together with the 6-jet reconstruct to a 



top quark. After such reconstruction, the signal events should show a mass peak over the 
WW background. Note that the WW "background" is itself the Higgs detection channel. 

Of course, ti background would always be present. The WW ti scattering via t- 
channel b quark exchange is suppressed in phase space compared to WW ^ h, H production 
followed by Higgs decay. When h, H ^ ti threshold opens up (not until 400 GeV or so), 
one would have genuine vvti ^ uu + bb + 4j background. These again can be distinguished 
from uutc production by event topology and jet counting. Since the ti/tc ratio is not that 
large [^H up to nih, h — 500 GeV, they do not pose a major threat. However, as seen from 
Figs. 2 and 3, for Higgs mass beyond 400-450 GeV or so, the signal cross section has also 
become too low and the WW background itself may start to become serious. 



C. Comparison of Different Processes 

It is of interest to point out the difference between vvtc production and other tc pro- 
duction processes. The e^e'^ ^ Z* —>■ tc, ic process, though rather clean, has very 
suppressed rate because the Z-t-c coupling is loop-induced (GIM mechanism is intact in 
the present model context). It is clear that e^e^ ^ Z* —>■ H{h)Z — > tcZ has identical 
sin^ a dependence as the WW fusion process. However, this production mechanism is less 
promising since it suffers from s-channel suppression (cross section decreasing as 1/s) at 
higher energies, and at the 500 GeV NLC, the rate is already a bit too low 0. 

The e"'"e~ Z* ^ h{H)A process is also s-channel suppressed, hence it is not particu- 
larly interesting at higher energies. But it does offer the intriguing signal of like-sign top 
quark pairs via h{H)A — > ttcc, ttcc, signaled by like-sign W plus bb events. Furthermore, 
the effects are the largest in this case when sin^ a ^ or 1, which is complementary to the 
sin^a domain of interest, Eq. (13), for the e+e^ i^e^^etc process. At the 500-600 GeV 
NLC, the rates for the two processes are comparable, both leading to only a handful of clean 
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events. Thus, though faUing short of making a definitive study, the 500-600 GeV NLC can 
cover the full range of sin^ a and offer us a glimpse of whether FCNH couplings exist or not. 

Turning away from e^e~ linear colliders, the process — > /i, -ff, A ^ tc [|18| at 
a possible future muon collider capitalizes on the larger Higgs-/i-/i coupling and a sharp 
Higgs resonance peak. However, because of the narrow width of the Higgs boson, this would 
demand precise tunings of the muon energies to find the Higgs resonance. In contrast, 
the beauty of the WW fusion process of Ref. M and discussed here is that no energy scan 



is necessary. It is not yet clear whether a high energy muon collider can be built or not [19 



Finally, let us compare with prospects at the LHC The challenge for VV — > /i, H ^ tc 
production search is the enormous background. It has been pointed out, however, that one 
might be able to directly probe for FCNH coupling strengths via the eg ^ tA ^ ttc, ttc 
production process at the LHC [0, which does not depend on sin a. Once again there is 
the intriguing signature of like-sign top quark pairs. The event rate is not very high since 
the raw cross section is at the 80 fb level |T^, and one still needs to make event selection 
cuts. Although promising, background rejection would certainly still be a major issue, as 
is almost always the case for interesting new physics at hadron colliders. In contrast to 
the high rate environment of the LHC, however, all high pt events at the NLC would be 
recorded and scrutinized. We stress that the search for FCNH effects via tc production is 
really part of the Higgs program. By studying the VV fusion processes alone, the relative 
large number of events in uutc mode (hundreds of events) and the concurrent study of 
Higgs boson properties via the W^W~ and ZZ modes (thousands and hundreds of events, 
respectively) should allow one to measure the h, H ^ tc branching ratios, which in turn 
can lead to a determination of the FCNH coupling. Thus, this has the advantage of being 
a complete program, and would be complementary to the eg tA process at the LHC. 
However, since it would only be fruitful for -^i > 1 TeV, the fulfillment of the program 
would certainly come after the studies at the LHC. 

At any rate, we expect the study of putc production via WW fusion to be quite feasible. 
We urge that a dedicated simulation study of this process for the NLC be carried out. 



V. CONCLUSION 

In summary, we have extended the work of Ref. [8] on tc-production via e~^e~ i^e^^etc 
at the NLC. We elucidate that the particularly promising mass range is when both and 
tuh are of order the weak scale. This is quite different from the parameter range discussed 
by the authors of Ref. [8] where one of the Higgs is taken to be as heavy as 1 TeV, and 
consequently the a^utc they obtained is smaller than ours. With Higgs masses in the range 
of 200-350 GeV, we find that a^utc could reach almost 10 fb. The sin^ a dependence is mild 
for 0.1 < sin^ a < 0.9, and a,^,ytc is greater than 1 fb as long as ^/s > 1 TeV. Given a 
significant cross section as such, this mode should be searched for carefully at future e^e~ 
Linear Colliders such as the NLC. 
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APPENDIX A: THE HELICITY METHOD 



The helicity method is particularly suited for numerical manipulations of scattering am- 
plitudes. For particles with spin, one constructs explicit representions for their helicity wave 
functions so that the relevant Feynman amplitudes can be written into numerical forms [jl2 . 
Consequently the squaring of scattering amplitudes may be performed numerically. 

For fermions, we choose the Weyl basis with the following represention of 7-matrices: 



7"= I,*!, f], (A2) 



or collectively 



.7^: 

with = (1, ^a). The chiral projection operator P± = (1 ± 75)/2 is then given by 

' °' p--(' ?). (A3) 



oy " vo 1 

where P+ and P_ project onto upper and lower components of Dirac four-spinors. 

In the Weyl basis, the Dirac spinor u{p, A) for a fermion with momentum p and helicity 
A is given by 



u{p, +) 



u4x = +)j 



«(p--).(::<^::))^(:;|;:|)^ <-> 

where u± = ±\p\ and \p±) denote the two-component eigenvectors of the helicity op- 
erator h = p- a/\p\ with 

cos ^ \ , „ , / — e"**^ sin ■ 



where 6 and are angles specifying the direction of p, i.e. 

p=\p\ (sin 9 cos 0, sin 9 sin 0, cos 9) . (A6) 
Similarly, the spinors of anti-fermions in the Weyl basis are given by 



v-iX = +)J \-^+\p-) 

We note that the helicity wave functions of spin 1 particles can be constructed out of the 
two building blocks: \p+) and \p—). 
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For any fermion(anti-fermion) line which contains arbitrary numbers of interaction ver- 
tices with bosons, the associated amphtude must be a hnear combination of the structures 

w^...Y-i''^^P-W2, (A8) 

and 

w^...YYl'"P+W2: (A9) 

where wi = i ) and W2 = \ I can be either u or v. Note that, for simphcity. 
Km- J \W2-J 

we do not specify the momentum and hehcity dependence of the spinors. Using the Weyl 
representations of Dirac spinors and 7-matrices, the above two structures are simphfied into 

«;I±---7^7-7>2-, (AlO) 

and 

w{^---Y_-fir-w2+, (All) 

where the sign in the subscript of depends on the number of 7-matrices inserted between 
the spinors. Using Eqs. (A4) and (A7), one can express Eqs. (AlO) and (All) as hnear 
combinations of 



(pi±|---7^7-7^|P2±), (A12) 

and 

{p,±\---Y-iir-\P2±). (A13) 

Defining the conjugate spinors as: 

|p±)=i(72(|p±))*, (A14) 

where 

\fr) = \f-) = 

= = (A15) 



Then, applying the relation 



'=F 

we have 



<y2{ilT<J2 = 72, (AI6) 



(pi ± I • • • 7^7-7^|P2±) = (P2 ± |7^^7;7^ • • • |Pi±), (A17) 

and 

(pi ± I • • • Y-lll'-\P2±) = (P2 ± |7^7-7^ • • • IPii)- (A18) 
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We now apply the above formalism to calculate e~^{pi)e~{p2) — * i^e{P3)^e{p4)t{pt)c{pc) . 
First, the amplitude for this process has been written in Eq. (4) with its fermionic part 
denoted as A - B ■ C. The explicit forms of A, B and C as shown in Eqs. (6) and (8) can be 
easily obtained by using Eqs. (A4)-(A18). Second, we note that the product A ■ B involves 
a contraction of Lorentz indices associated with matrices 7+ and 7^_. Such contractions can 
be evaluated easily via the "Fierz-like" relation 

{1+)ij{1^,-)kl = {l-)ij{'l^i+)ki = '2,5ii5kj (A19) 

with i ,j, k and / being indices in spinor space. Indeed, from (A19), we have 

{P3 + Mpi+) ■ {Pa - \l^+\p2-) = 2{p3 + \p2-) {pa - \pi+) (A20) 

From Eqs. (4), (6) and (8) and Eq. (A20), one now has the full helicity amplitudes iM{Xt, Ac) 
for e+(pi)e~(p2) — ^ i^e{P3)t^e{p4)t{pt)c{pc) , which can be easily incorporated into the numer- 
ical program ONETOP 



APPENDIX B: THE INTERFERENCE OF FEYNMAN AMPLITUDES 

In this Appendix, we discuss the interference effects of scattering amplitudes arising from 
different neutral Higgs bosons. Let us use iA4^ to denote the amplitudes of e+(pi)e~(p2) — *■ 
^e{P3)^e{P4)t{pt)c{pc) coutributious from neutral Higgs bosons S = H and h. That is 

zA^^ = zM(e+e- Uei^eS*{q)) x — x iM{S*(q) ^ tc), (Bl) 

— + imsTs 

where S*{q) denotes the off-shell S with momentum q. The total cross section (Tuutc = 
cr(e+(pi)e"(p2) h'e{P3)i^e{PA)t{pt)c{Pc)), IS giveu by 

{2'nY f d^P3 d?PA d^Pt d^Pc 



2s J (27r)32E3 {27t)^2E4 {2T:f2Et {2'Kf2Ec 
X 5\pi+p2-p3-PA-Pt-Pc)\iM'' + iM''\\ (B2) 

One can separate (Tuvtc into diagonal and interference terms, i.e. 

In the narrow width limit H{h) ^H{h)-, it is well known that 

= cT(e+e- VeVeH{h)) x BR{H{h) tc). (B4) 
However the interference term a^^^ is more complicated. From Eqs. (B1)-(B3), we obtain 
.H-h _ „„„2 „, „,„2 „ (27r)^ [ d^ps d^p4 d^q f 3 ^ 1 



0",/,/^;= = — COS a sin a / - — — — — — / (2tt) dq 

""^^ 2s J {2n)^2E3 (27r)32E4 (27r)^ J ^ ' 2Eq 

X +P2-P3-PA-q) \iM{e'^e' -> Ve^eHsyCiq)? 

^ 2Re [- — — ) 

\q^ - mjj + zttihTh q^ - ml- irUhTh ) 
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where Eq = \/q'^ + q-^, and HsM*{q) has been used to replace H*{q) or h*{q) since we have 
factored out the mixing-angle dependence cos^ a sin^ a. 

In general, H and h are not degenerate. Without loss of generality, we may assume 
niH > ruh so that m% — m\ — L x mhTh with L > 0. Furthermore let us take x = q^ — nijj. 
With a little algebra, the propagator part of Eq. (B5) can be written as 

1 1 



J{2nfdq'^---2Re 
^Jdx^ 

-I 



rrijj + irriH^H q^ — rni 




27r \x^ + mlrl {x + LruhThy + mlrl 
, 1 jmlrl - 2mHrHm,r, + (L^ + l)m^r^) , 



If H and h are precisely degenerate, i.e. L = 0, and sin^ « = 1/2 which implies = T^, 
the second term on the r.h.s. of Eq. (B6) vanishes, while the first term eventually gives rise 
to —{(T^ytc + ^vvt^ which cancells completely the diagonal contributions as expected. 

For L > 3, we note that the (L^ + 1) m^F^ term in Eq. (B6) already dominates over both 
m^Vlj and mfj'\, provided F^ is not overly suppressed by too small a sin^ a. Therefore, in 
Eq. (B6), one may neglect the combination m^T^ — 2mi{V HiiihTh with respect to (L^ + 
l)m^F^. In this approximation, one can show that the two terms on the r.h.s of Eq. (B6) 
lead to a vanishing interference term cr^^j* in the narrow width limit of ^H{h) ~^ 0. 

To see this, note that in the limit of '^h(}i) 0, the dominant contributions to the 
X integration in Eq. (B6) comes from the vicinities of x = —LrrihTh and x = 0. For 
X ~ —LmhTh, we discard the term l/(x^ + m|^F|^) on the r.h.s. of (B6) while the remaining 
terms are rearranged as follows: 
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/ 
/ 
/ 



1 1 

"^"^2^ ■ ■ ■ {x + LmHT^y + mlTl^'^''> 

''''271 {x' + m%rl) {{x + Lm,F,)2 + m^F^) ^ 

7r{x + Lm,r,r + mlrl2m,r,^^''^ ' 

, 1 {L' + 1) . m,F, ^r2^^4 .R7^ 

''^■■■7r(x2 + ml,Fl,)((x + Lm,F,)2 + m|F2)2m,F/ ' ^ ' 



In the limit F/j — > 0, the first term on the r.h.s. of (B7) can be simphfied by 

1 mhTh 



n {x + LmhVhf + mlTl 



5{x + LmhVh), (B8) 



while the factor l/2mhTh leads to BR(/i — > tc) when combined with other terms in Eq. 
(B5). Similarly, the second term on the r.h.s. of (B7) can be simplified by 

1 (L^ + 1) mlVl . mnrh ^ (L^ + l)m|F| 

TT {x^ + m%rj,) {{x + Lm,F,)^ + m|F^) {L^min + mlT%) '^"^ + ^'^ 

~ Six + LmhTh), (B9) 



15 



where we have used ^ 1. Clearly the two terms on the r.h.s. of (B7) cancel completely, 
hence (T^^e* receives no contribution from x —LrrihTh. By similar arguments the integra- 
tion region a; ~ also gives no contributions to cr^ytz- We therefore conclude that a^^^ = 
provided we neglect the combination mj^Tj^ — Hf^hXh with respect to (L^ + l)m\T\. 

Keeping the myr'\j — 2mH^HTnhrh term in Eq. (B6), the resulting cr^^c'* is at least 0(1/ L'^) 
suppressed compared to the total diagonal cross section ((7;^^^ + crl^,^tc)- 
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FIGURES 



FIG. 1. Feynman diagram contributing to the process e^e — > Vf,Vf,tc. 

FIG. 2. The cross section Oyytc as a function of nih with vhh = 1 TeV and sin^ a = 1/2 for 
= 0.5, 1, 1.5, 2 TeV (bottom to top). SoUd Unes are for the fuh calculation, while dashed 
lines are from Ref. [8] which uses the effective W approximation. 

FIG. 3. The cross section Ovvtc as a function of m/j with mu = 300 GeV and sin^ a = 1/2 for 
sV2 = 0.5, 1, 1.5, 2 TeV (bottom to top). 

FIG. 4. The cross setion auuHsM of e'^e~ — Ve^eHsu as a function of mn^m 
sV2 = 0.5, 1, 1.5, 2 TeV (bottom to top). 

FIG. 5. The branching ratio BR(/i — > tc) as a function of mh for sin^ a = 0.1, 0.3, 0.5, 0.7, 0.9 
(top to bottom). 

FIG. 6. The effective fraction cos^ aBK{H —>■ tc) + sin^ aBR(/i tc) as a function of sin^ a 
with mH = 250 GeV and nih = 240 GeV. 

FIG. 7. The cross section a^i^tc as a function of sin^ a with mn = 350 GeV and ruh = 200 GeV 
for sV2 = 0.5, 1, 1.5, 2 TeV (bottom to top). 
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